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The recently reported ultracold trapped Fermi gases with tunable atomic scattering 
length [1-19] in the vicinity of a Feshbach resonance stimulated a large number of theoretical 
investigations. In this letter the dynamics of strongly interacting dilute Fermi gases (dilute 
in the sense that the range of interatomic potential is small compared with inter-particle 
spacing ) is investigated in the framework of the hydrodynamic theory [20-29]. 

Instead of using the scaling approximation of Refs[20-22, 25,27], the polytropic approx- 
imation of Refs. [24,28] or the perturbative approximation of Ref. [26], we will construct a 
rigorous upper bounds to the collective frequencies. 

Our starting point is the variational formulation of the Kohn-Sham time-dependent the- 
ory for a dilute two-component Fermi gas in a trap potential V ext (r) = (m/2)(uj 2 _(x 2 + y 2 ) + 

5 J dt <V\ihd t - H\V >=0, (1) 

where |\P > is a product of two Slater determinants, one for each internal state built up by 
the Kohn-Sham orbitals ipi, and H = T+U is Hamiltonian in the local density approximation 
(LDA). 

We consider two approximations: 

(i) local transform ipi 0jexp(ix), where ipi and x are re& l functions, 

and 

(ii) < <j>\T\<j> >^ (l/N)J(t TF (n)+t w (n))n(f,t)d 3 r, 

where n is the density, normalized to the total number of atoms N, J nd 3 r = N, [0 > is the 
product of two Slater determinants built on 0, alone, ^tf(^) = (37r 2 ) 2 / 3 (3/i 2 /(10m))n(r, t) 2 ^ 
is the Thomas-Fermi (TF) kinetic energy density, and t w (n) = (fr 2 /(8m))(Vn) 2 /n is the von 
Weizsacker kinetic energy density. 

Using the approximations (i) and (ii), we obtain from Eq.(l) [20-22] 

dib h 2 

where ip = y^nexp(ix), V xc (f,t) = [ dn ^ ] n =n(f,t) an d e(n) is the ground state energy per 
particle of a uniform system. The only difference from equations holding for bosons [30,31] 
is given by density dependence of e(n). 
It is useful to rewrite Eq. (2) in a form 

% + V(nv) = 0, (3) 



dv 1 „ /T , d(ne(n)) 1 9 ft 1 ~ „ , 
« + m V(,/ - + "4^ + 2 mV - 2^ VH V ^ = °' < 4) 
where v = (K/m)Vx is the velocity field. 

It was shown in Refs. [20,21] that for experimental conditions of Refs. [14,16,17] the quan- 
tum pressure term in Eq.(4) can be neglected. For the reminder of this Letter we will use 
this hydrodynamic approximation. 

For the negative S-wave scattering length between the two fermionic species, a < 0, in 
the low-density regime, kp \ a \<ti 1, the ground state energy per particle , e(n), is well 
represented by an expansion in power of kp \ a | [32] 

e(n) = 2E F [^- - -^—kp | a | +0.055661(A; F | a |) 2 - 0.00914(A; F | a |) 3 + ...], (5) 
10 on 

where E F = h 2 kp/(2m) and kp = (S^n) 1 ^ 3 . In the opposite regime, a — > — oo (the Bertsch 
many-body problem, quoted in Refs. [33]), e(n) is proportional to that of the non-interacting 
Fermi gas 

where the universal parameter (3 [10] is estimated to be (3 = —0.56 [34]. The universal limit 
[10,34-38] is valid at least in the case where the width of the Feshbach resonance is large 
compared to the Fermi energy as in the cases of 6 Li and m K. 

In the a — > +0 limit the system reduces to the dilute Bose gas of dimers 



e(n) = E F (-l/(k F a) 2 + a m k F /(Qn) + ...,) (7) 

where a m is the boson-boson scattering length, a m w 0.6a [40]. 

A simple interpolation of the form e(n) ~ EpP(kpa) with a smooth function P(x) was 
considered in several papers. In Ref.[20] it has been proposed a [2/2] Pade approximant for 
the function P(x) for the the case of negative a 

P(x) = --2 - ' ,' , (8) 



5 

where S 1 = 0.106103, S 2 = 0.187515, 5 3 = 2.29188, 5 A = 1.11616. Eq.(8) is constructed 
to reproduce the first four terms of the expansion (5) in the low-density regime and also 
to exactly reproduce results of the recent Monte Carlo calculations [34], (3 = —0.56, in the 
unitary limit, kpa — > — oo. 
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For the positive a case ( the interaction is strong enough to form bound molecules with 
energy E mol ) we have considered in Ref.[21] a [2/2] Pade approximant 

P( x ) = E^L + aiX ± a2%2 (g) 
1 J 2E F l + a 3 x + a 4 x 2 ' 1 ' 

where parameters a are fixed by two continuity conditions at large x, 1/x — > 0, and by 
two continuity conditions at small x, ax = 0.0316621, a 2 = 0.0111816, a 3 = 0.200149, and 
a 4 = 0.0423545. 

In Ref.[41] a Pade approximation has been considered for the chemical potential. Authors 
of Ref. [28] have used a model for P(x) , interpolating the Monte Carlo results of Ref . [36] across 
the unitary limit and limiting behaviors for small \x\. We note here also the BCS mean-field 
calculations of Ref. [27] 

The hydrodynamic equations after linearization take the form 

w Sn+ ^ (nV( ^ Jl6n)) = ' (10) 

where 5n(r, t) is the change in the density profile with respect to the equilibrium configura- 
tion. If we consider oscillations with time dependence 5n oc exp(icut), Eq.(10) can be reduced 
to a Hermitian equation [26] 

1 (Hi 

an 1 m 

where / = d Sn, and the equilibrium density, n, is given by equation 

, = V eIt+ d -^i, (12) 

where /i is the chemical potential, in the region where n(f) is positive and n(r) = outside 
this region. 

The Hilleraas-Undheim method [42] for upper bounds to eigenenergies of excited states 
of atoms can be generalized to the case of Eq.(ll). Defining the functional I[x] by 

I[ X ] = -(1/m) J X V{nVx)d"r/ J x f^^ Y^r, 

we get /[/] = uj 2 , where / is the solution of Eq.(ll). 

Let the sequence of eigenvalues of Eq. (11) be denoted by uo\ < uo\ < ...and let Xi>X2,---,X ? 
form a set of q linear independent functions. Introducing the function x — J2l=i c rXr, where 
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Ci, c 2 , ...c g are g variable parameters, we see that the functional I is stationary if dl/dc r = 0. 
Then the values of the q roots, Co 2 < G)\ < Co 2 < ... < Co 2 , of the secular equations 




provide upper bounds to the q lowest eigenvalues of Eq.(ll) u 2 > uj 2 for r = 1,2, ...,q. It 
follows that as extra terms are added to the expansion of Xi a given root of the secular 
equation u r is decreased. 

For the most interested case of M = modes, we can put in Eq.(13) q = 3, Xi — 1> 
X2 = (^ 2 + y 2 ) and x 3 = which give 



where rj = (3 + 4A 2 )C - (3 + 6A 2 ), C = W^f, J « = fx l n(x)dx, x = ^x 2 + y 2 + X 2 z 2 , 



A = l} z /u]_, and ± signs refer to the transverse and axial mode, respectively. 
To calculate (, we have used the following expansion 

n(f) w (1 - /JV ext (f)) 1/(2 - rt Eci^U* 5 )]*, (15) 

i=0 

where parameters (3, p and q are fixed by requiring that n(f) must satisfy a variational 
principle 5 J n(V ext + e(n))d 3 r = with a subsidiary condition / nd 3 r = N. 

From the Table 1 one can see a very fast convergence for ( and consequently for u±, the 
first term of the expansion (15) has accuracy 1CT 3 [43], while the sum of the first three terms 
has accuracy 1CT 4 . 

It is easy to show that our upper bounds, Eq. (14), give exact solutions for frequencies 
of the breathing modes for the polytropic equation of state, e(n) ~ n 1 . 

The hydrodynamic equation is expected to be applicable for describing the macroscopic 
excitations of the system up to energies of the order of the energy gap, A, needed to break- 
up a Cooper pair. But for the trapped gas, A is a function of r (A decreasing when we 
go away from the center). It is natural to assume that the condition of the applicability of 
hydrodynamics to describe the macroscopic excitations of the system at T = is [22,44] 

Ml , . 

x « 1, (16) 

where fl is the frequency of the macroscopic excitations and the mean energy gap is given by 
A = / n(r)A(r)df/N. To calculate A we have used results of Refs. [45,46]. The predictions of 




r ? ± V / 7 ? 2 -8A 2 C(5C-9) 
(5C - 9) 



(14) 
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Eq.(14) for the radial breathing mode frequency of the cloud of 6 Li atoms are shown in Fig.l. 
For e(n) we have used two approximations, the Pade [2/2] approximation of Refs. [20,21] and 
the parameterization of Ref.[28]. It can be seen a very small difference between these two 
appr oximat ions . 

In Fig.l, we have also compared the hydrodynamic predictions with experimental data 
[17]. We have used Ref. [47] to convert a to the magnetic field B. Near the unitary limit there 
is a very good agreement with experimental data [17]. We note here that two experimental 
results [14,17] and [16] for oo + are still about 10% in disagreement with each other, which 
is not fully understood yet. At [iV^a/a^J -1 = —1.26 the hydrodynamic upper bound 
begins to be below the sharp increased measured frequency that confirms a breakdown of 
hydrodynamic theory [17]. 

In Fig. 2, the calculated hydrodynamic upper bounds to the axial breathing mode fre- 
quency of the cloud of 6 Li is compared with experimental data [16] in the BCS-BEC crossover 
region. It can be seen that the difference between the two approximations for e(n) is prac- 
tically negligible and both approximations give a very good agreement between calculations 
and experimental data [16] for [^^a/aho]' 1 > -1.14. However at [N^a/aho]- 1 = -1.31 
the hydrodynamic upper bound begins to be below the measured frequency that indicates, 
as in the radial mode case, breakdown of hydrodynamics, even though for the axial mode 
case the ratio between the collective energy and the gap energy, Eq.(16) is very small, 
Tiuj-/K ps 0.01, at [N^-^a/aho]^ 1 = —1.31 [48]. This breakdown of hydrodynamic theory 
may be related to the finite-temperature effects. 

In summary, we have constructed the rigorous upper bounds to the collective frequencies 
of the Fermi gas under harmonic confinement in the framework of the hydrodynamic theory. 
The bounds are compared to experimental data on confined vapors of 6 Li atoms. It is 
shown that the hydrodynamic upper bound begins to be below the measured frequency at 
[A^ 1 / 6 a/a/ lo ] _1 Ps —1.3, although the ratio between the collective energy and the gap energy, 
Eq.(16), for the axial mode is very small. 

We thank J. Thomas for stimulating this work and for sharing with us the experimental 
results prior to publication. We also thank R. Grim and M. Bartenstein for providing us 
with the updated experimental data. 
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Table 1. ( in the BCS region as a function of the dimensional parameter X = (A rl / 6 a/a ?i0 )" 
and the number of terms, /, in the expansion (15). The [2/2] Pade approximation of 
Refs. [20,21] is used for the energy per particle e(n). 



-0.1 2.2587 2.2572 2.2577 2.2575 2.2576 

-0.2 2.2658 2.2637 2.2642 2.2640 2.2641 

-0.4 2.2753 2.2732 2.2735 2.2735 2.2735 

-0.6 2.2801 2.2789 2.2787 2.2789 2.2788 

-0.8 2.2821 2.2818 2.2814 2.2816 2.2815 

-1.1 2.2822 2.2830 2.2825 2.2827 2.2826 

-1.4 2.2810 2.2825 2.2819 2.2821 2.2820 

-1.7 2.2792 2.2811 2.2805 2.2807 2.2806 

-2.0 2.2774 2.2795 2.2789 2.2791 2.2790 

-3.0 2.2720 2.2742 2.2738 2.2739 2.2739 
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Fig. 1. Radial breathing mode frequency ui + in units of u± as a function of the di- 
mensional parameter X = (A^ 1//6 a/a/j ) _1 . In the unitary limit, a — > — oo, one expect 
u/u± = y^lO/3 ~ 1.826. The solid line and the dashed line represent the hydrodynamical 
upper bounds calculated using for the energy per particle the [2/2] Pade approximation of 
Refs. [20,21] and the parameterization of Ref.[28], respectively. The circular dots with error 
bars are the experimental results given by the Duke University group [17]. 
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Fig. 2. Axial breathing mode frequency u- in units of uo z as a function of the dimensional 
parameter X = (A^ 1//6 a/a/ lo ) _1 . The circular dots and diamonds represent the hydrodynam- 
ical upper bounds calculated using for the energy per particle the [2/2] Pade approximation 
of Refs. [20,21] and the parameterization of Ref. [28], respectively. The triangles with error 
bars are the experimental results given by the Innsbruck group [16]. 
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